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EQUATIONS OF MOTION FOR THE X-14 AIRCRAFT 


Introduction | 

l 

The research described in this report presents partial results of 
a study of the control and power requirements of the X-14 VTOL aircraft 
under NASA Grant NGR-05-004-051 . 

The complete equations of motion for the X-14 are derived in this 
report. The fundamental assumption is that the aircraft is a single 
rigid body. The equations of motion are derived x^ith respect to a set 
of axes fixed to the aircraft. Additional assumptions used are that any 
wind disturbances are irrotational that the tx^in engines used on the air- 
craft rotate in the same direction at the same speed and that the engine 
exhaust is diverted by means of vanes to provide a direction varying thrust 
vector. 

The equations obtained are subsequently linearized about various re- 
ference conditions and numerical values for the trim parameters and the 
stability derivatives at these conditions are tabulated at the end of the 


report . 
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CHAPTER I 

DEVE L OPMENT OF THE EQU ATIONS OF MOTION FOR THE X-14 A IRC FAFT 

A. General Equations of Motion for VTOL Aircraft 

This development of the equations of motion for a! VTOL aircraft 

i 

assumes the aircraft is rigid and the origin of the coordinate system 
is at the mass center [1, 2, 3]. Additional assumptions used in this 
development of the equations of motion are: 

1. The VTOL aircraft has diverters and vanes affecting the engine 
exhaus t 

2. Any wind disturbance is irrotational 

3. Engines rotate in the same direction at the same speed. 

A standard aircraft body axis system is used where the aircraft is 
symmetric about the x-z plane with the positive y-axis pointing cut the 
right wing, the z-axis down, and the x-axis in the forward direction of 
flight (stability axes). This is a right-handed coordinate system. 


For 

the system, 

the velocities, 

angles 

, angular velocities, force 

and moments associa 

ted with each axis are 

shown in Table 

A- 1-1 and 

Figure 

A- 1-1 [4] . 

Table A- 1-1 




Axis 

Linear 

Velocity 

Force 

Angle 

Angular 

Velocity 

Moment 

X 

u 

X or F 

X 

<}> 

p 

L or l 

y 

V 

Y or F 

y 

9 

q 

M or m 

z 

V 

Z or F 

z 

Y 

r 

N or n 
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From Etkin 11 J, the time derivative of the angular quantities can 
be expressed as 

« 

0 = q cos r sin 

<£ *= p + q sin $ tan 0 + r cos £ tan 0 j 

T = (q sin $ + r cos $ ) / cos 9 j 

Since the analysis is for a VTOL aircraft, it becomes necessary to 
introduce some notation not normally found in airplane equations. A 
listing of this additional notation is in Figure A- 1-2. 


Inertial Force s 

In writing the equations of motion, all inertial components are 

expressed as D'Alembert forces, i.e., 

EF . + EF. . « 0 

external inertial 


where 


EF = -ma. 

inertial 

Using this notation the inertial forces along the three axes are 


xi 


-m (u + wq - vr) 


F . = -m (v - wp + ur) 
yi 


zi 


-in (w - uq + vp). 


Likewise, the moments are 


M 


xi 


M 


yi 


M 


zi 


_ Pl + qr (I - I ) + (r + pq) I 

x y z xz 

- q I + pr (I - I ) + (r 2 - p 2 ) I 

y z x r xz 

- r I + pq (I - I ) + (p - qr) I 

z r x y 1 H xz 
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Figure A-l-2 

Definitions of Non-standard Notation 

1, w subscript is for wind velocity 

2 , «* jet exhaust velocity from the diverters 

3 # o m thrust diverter angle measured in the vcrticle plane from the 
positive z-axis 

4. X sc exhaust side vane angle measured positive about the x-axis by 
the right hand rule 

5. U n - [(u - u ) + (w - v ) ] 

U w w 

— 2 , v 2 , >2, 1/2 

6. U Q - [(u - u w ) + (w - w w ) + (v - v w ) ] 

7. a = angle of attack; sin ° *= (w - w u )/ u 0 

8. 8 “ side slip angle; sin 3 ■ (v - v^J/Uq 

9. ** distance engine intake is in front of tne center of gravity 
19. z = distance engine intake is below the center of gravity 

11. x 2 = distance engine exhaust pivot point is in front of the center of 
gravity 

12. z 2 “ distance engine exhaust pivot point is below the center of gravity 

13. is the effective length where the exhaust will impinge on the diverter 
vane from the end of the jet engine 

14. I *= moment of inertia of one engine about the axis of rotation 

e 

15. fi ■ angular velocity of an engine measured positive in the direction 

e 

of the x-axis 

16. N.E. *» number of engines 

17. T - total thrust from the diverters 

18. L M , H apparent aerodynamic damping near hover not directly 

p q r ‘ 

attributable to the effects caused by C £ * c m * c n 

p q r 

19. Lr , M , N Moment effects due to reaction control nozzles 

^ O A 

x y °z 


*=j et velocity at the engine exhaust 

total installed jet thrust 
longitudinal diverter efficiency factor 
lateral diverter vane efficiency factor 
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where 


I & /( y 2 + Z 2 )dm, I « /(x 2 + z 2 )dm, I z - /(x 2 + y 2 )dtn 

^ ^ A 

I £ /xzdm I - /xydn I - /y^m 

XZ / 


with 

X “0 and I “0 

xy yz 

8lnce the aircraft is symmetric .about the xz plane. 

i 

In addition to aircraft inertial terms, there is a contribution from 
the rotation of the engines to the moment equations. If the engine speed 
is almost constant, the contribution to the moment equations can be expressed 

as 

M « 0 

e 

x 

H « (N.E.) I n r 

e e e 

y 

m - -(n.e.) i n q 

e e e 

z 

Gravitational F o rces 

Since the origin of the coordinate system is at the center of gravity, 
the gravitational forces do not enter into the moment equations. Their 
contribution to the force equations is 


-W sin 9 


g 


F ** W sin $ cos 6 

y g 

F «* W cos <J> cos 6 
z 

g 


Mass Flow Effects 


There are two contributions to mass flow effects, one at the engine 
intakes and a second effect at the engine exhausts. If the mass of the fuel 
burned is neglected, the forces due to the mass flow can be expressed as 


F 

X. 

m 

» -m 

(u - 

u 

w 

b 

sin 

a cos 

F 

y m 

■ -m 

(v - 

v - 
w 

b 

sin 

*> 

F z. 

m 

*» -in 

(w - 

V + 

w 

b 

cos 

0 cos 
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where iri is the entering mass flow rate which is approximately equal to the 
mass flow rate out, is the relative exhaust velocity of the jet, ° is 
the diverter angle, and X is the side vane angle. Recalling that the gross 
thrust equals mass flow rate times relative velocity, 

T =■ m Vj 


or 

T 

Upon making this substitution, the mass flow contributions to the equations 
of motion can be written as 

F E T [sin a cos X - (u - u )/V ] 

x. w j 

m 

F *= T [sin X - (v - v )/ V ] 

y. w j 

in 


m 


=-T [cos a cos X + (w - w )/V ] 

w j 


The mass flow contributions to the moment equations enter in because the 
intakes and exhausts are not at the center of gravity. The moments are 


M 


x. 

m 


= T z (v - v )/V 

X W j 


in 


M 


xa 


T [x x (w - v w ) " z 1 (u - u w )]/V 


in 


M 


ra 


'"T x (v - v )/V 
i w j 


in 


where x^ and z^ are the distances from the center of gravity to the center 
of the intakes. The moments caused by the exhaust are 


M 


m 


out 


M 


m 


cut 


M 


m 


= -T sin X(~2 + cos cr) 


T cos X (z sin a + x^ cos °) 
T sin X (x„ - 2 .^ sin a) 


out 


where X2 and Z2 are the distances from the center of gravity to the center 
of the exhausts, and 2 ^ is the effective length where the exhaust will 


impinge on the diverter vane. 
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The total moments caused by engine mass flow terms can now be 
expressed as 

M ■ T f-sin X (z. + l. cos o) + z (v - v )/V ] 
x. L 1 X w j 


m 


M 


T { cos X (z 2 sin o + * 2 cos a) + (w - w^) - z^ Cxi— u^) ] /V } 


in 


M z *= T [sin X(x 2 - ^ sin c) - (v - v w )/V ) 
m 

Additional moments are caused by discharge nozzles in the tail and on the 

wing tips. These reaction control moments can be written as 

M 
x 

M 

y 

M 


«= L r <5 

X R.C. 5 * x 


H- 6 

y R.C. °y y 

N. <5 
o z 


R.C. 

where <5 , 5 and 5 are the reaction nozzle openings expressed in degrees, 
x y z 

Aerod y namj c Ex pre ssi ons 

The complete aerodynamic expressions for a conventional airplane can 
be written as 

x ■= f u 2 s c 

aero ^ o x 

Y = •£ U 2 S C 
aero 2 o y 

z = ~ u 2 s c 

aero 2 o z 

■£ U 2 S b c„ 


aero 2 o 


M 


aero 2 o 


-£ U 2 S c C 


m 


N 


* S u 2 S b C 


aero 2 o 

where C . C , C , Co , C and C are the non-dimensional aerodynamic 
x y z m n 

r , .2 , . . 2 . 1/2 

coefficients, p is the air density, U q » l(u - u^) + (w - w^J J , 

U - [(u - u ) 2 + (w - w ) 2 + (v - v ) 2 } 1 ^ 2 , S is the wing area, b is the 
o w w w 

span and c is the wing chord. 
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Now use a Taylor series expansion of the non-dimensional aerodynamic 

forces retaining only the first terms in the series [1, 2] and let the 

A 

vehicle velocities relative to the air mass be denoted by u ^ ^ u - u^, 

A A 

v . *» v - v and w , * w - w . Then, 
rel w rel w 


3C 9C 


3C 


C ** C + a + irrAl + 8 

x x-, rel e 

0 c 


BC 

c + — £ 


3C 


9C 


9C 


8 + P + ~ r + — - 6 

93 op ‘or 96 r 


9C 


9C 


C = C + •— ~ a + u 
z Zq 9a 


9C 


9C 


9C 


?. 7 . • Z - 

u + “ — q + ~~ — a + "vr d 
9u rel 9q 4 o a 06 e 


9C, 


sc ac ac 4 sc t 

c. - c. + — 6 + p + — r + - 53 - « + 3 — e T 

K a r 


'l ' 38 


9C 9C 

C = C + — 1 a + -t- 2 u 
in Dq 9 a 


9C 


9C 


m , m • , m - 

j ~ — q + a 9- — — 6 

ou rel 3q n oa 9o e e 


9C 


9C 


+ 

3C 


ac 


sc 


C n = C n + 33 B + 9p p + 9r r + 36 5 a + 36 5 r 
u “0 K a r 


where a is the angle of attack (sin a = w / U Q ) , 3 is the side slip 
angle (cos 3 = U Q / U ) , and 6 £ , 6 g and 6^ are the elevator, aileron and 
rudder deflections, respectively. 

Following Etkin [2], if the partial derivatives are written as 
non-dimensional stability derivatives, the aerodynamic portions of the 
equations of motion become 

X x^u 2 S[C +C ct+^C u ,+ C 5 1 

aero 2 o 1 x x u U x rel x e 
o a o u 6^ 

e 


aero 


J U 2 S|C +C 3+-C P + -C r+C 4] 

20 y o y 3 2tf y p 2U y r y 4 r 

o r o r 


aero 


^J. 2 S[C. + C. a + i C z u rel +2 § C , + 2 § C z _4+C z 5 C) 


2 o ‘ z z 

o a o u 


o q 


o a 5 


\ 
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aero 


fU o 2 SMC +C S + *C P + = C t r + C *. + 

2 0 l o 2U p 2 U q 4 r 5 fl 

c „ « ] ° ° 

*6 r 


M 


aero 


“ U 2 Sc[C + C a+ h C u + ^ c q + 
2 o m in U m rel/u m 


o u 


o q 


J- Q a-f C 6 ] 
2U V m 6 e 

e 


N 


aero 


r U 2 Sb[C + C 3 + - ' C p + ~ C r + C 
2 o n n 7 oT 7 n _ oTT n ,- n , 


2U “p 2U r 

o o a 


+ C <5 ] 


nr 


In addition to those stability derivatives normally found in airplane 
equations, for VIOL aircraft it is desirable to add the rotational 
damping observed on a hovering airplane. This can be expressed as: 

L aero VIOL L p P 

^aero VIOL ‘*q ^ 

N aero VTOL " N r r 

After sunning the component parts, the complete non-linear equations of 
motion for a VTOL aircraft can be written as in Fig. A-l-3. 
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Figure A-l-3 

Summary of Complete Non-linear Equations of Motion for a VTOL Aircraft 


T 1 

u « vr - wq - g sin 0 + — [sin o cos X - — (u - u )] 

^ c m V j W 

+ — pU 2 S [C +C ct + -jp C (u - u ) + C 6 ] 

2ra 0 x- x U. x w' x. e 

0 a 0 u o 

e - 2 

T 1 p ^0 ^ 

v ® wp - ur + g sin $ cos 0 1 — [sin X - (v -v ) ] + 

m v . w 


j 

[C Y + c s+ ^r-c P + V- C Y r+C Y M 
0 *3 2U 0 x p 2U 0 r x 6 r 


2m 


(A-l-1) 


(A-l-2) 


T 1 

w * uq - vp + g cos 6 cos 9 - — [cos a cos X + — — (w - w )] 
1 r " ' m V ■ w 


pU c l 

+ — ^ — - [C + C a + ~ C (u - u ) + v— C q 

2m 1 z n z b n z w 2U z 

0 a 0 u 0 q 


(A-l-3) 


+ ■—*- C a+ C o ] 
21' z ■ z. e 

0 a o 

e 


(I - 1 ) I _ 

p = - qr + ~~ (r + pq) + y [-sin * ( z 2 + l ± cos a ^ 

X X x 


Z, L h 0 2 Sb 

+ v: <v • V + i* > + r* { * + 0 n — lc i 0 + % 8 

j XXX P 


(A- 1-4) 


+ ^- C 


2U 0 p 


b 

2U 


1 P+fr-C, r + C 5 a + c 5 r I 


0 


(I - 1 ) 


q ■ 


z x , xz , 2 2, 1 ^ . T 


pr + — — (r - p ) - 7 - (N.E.) I « r + 


e e 


{ cos X (z. sin o + x„ cos a ) + — [x. (w - w ) - z. (u - u ) ]} 
l Z Vj l w i w 
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M V |JV~ <i 

+ -f q + — { X + — f 1C„ + C. a + ±-C„ <u - uj 


H, 


pu o 2sc 


I 


21 


,n 0 in a 


U n m 
0 u 


w 


(A- 1-5) 


+ T7T- c q + §7T- c a + c 5 ] 

2U 0 m q 2U 0 "6 m 6 6 


r = 


— y ~~ pq + l~~ ^ ~ qt) + i” I e n e q 

z z z 


N 


+ ~ [sin X (x - * sin a) - ~ (v - v ) ] + ~ r 

z j z 


(A-l-6) 


N r r, 2 , 

<5 p U Sb , i 

+ —5. 6 + — — [C + C 6 + C p + — C r 

l z * 21 z "0 "S 2U„ "p 2U Q n r 


+ C 6 + C 6 

n 6 2 n 6 

a r 
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B. Application of the Equations of Motion to the X~14 Aircraft 

Several of the above aerodynamic stability derivatives have a very 
small effect in the representation of the dynamics of the X-14. The 
determination of which stability derivatives have this insignificant 
effect was based upon examination of wind tunnel data and a knowledge 
of the physical airplane. The final justification for retaining or 
neglecting different stability derivatives comes from the fact that an 
analytical model was obtained which when used in simulator studies closely 
approximated the physical airplane [5]. Using these findings the above 
aerodynamic contribution to the equations of motion can be expressed as 


X » U 2 S (C T sin a - c n cos a) 

aero 2 o L D 

Y - f- U 2 S C 
aero 2 o Y 

2 

Z « ~£U S (C_ cos a + C sin a) 
aero 2 0 ^ ^ 

L “ | U 2 Sb (C. + C. 6 ) + L p 

aero 2 o * a p 

o 0 

a 

M ' - % U 2 Sc (C + C a+ C q + C 6 ) + M q 

aero 2 o n m 2U m m- e q 

o ct 0 0 o 

n e 

N U 2 Sb (C + C 6+C 6 ) + N r 

aero 2 o n n f r n. a r 

° 6 r 6 a 


where is the non-dimensional lift coefficient and is the drag 
coefficient [2] . 


Continuing to follow the development in reference [5], some of the 
above "power on" stability derivatives are expressed as "power off" plus 
11 A 11 terms where the 11 A 11 represents the difference between power on and 
power off effects. The aerodynamic portion of the equations of motion 
now becomes 

X - T U 2 S[(C, + & C. ) sin Ct - C cos a] 

aero 2 o L , „ ^ ce 1* u 


Y » ~ U I 2 S C v 
aero 2 o Y 


aero 


- tt U 2 S [(C T + AC,) cos a + C sin a] 

2 O L r r 

power off 


L = ~ U 2 Sb (C + C 6 ) + L p 
aero 2 o 4 4. a p 

o o 

a 


M 


aero 


%■ u 2 sc [ (c 

n 


m 


+ AC,.) + (C 

i-'l It 


power off 


+ A c ) ct 

\x ^ 

power off 


^ 7 - C q + C 0 ] + M q 
2U m iur e q 

o q 0 


N = | U 2 Sb (C + C 5 + C 6 ) + N r 

aero 2o n nr n_ a r 

o 6„ 


Upon incorporating these expressions for the aerodynamic forces in the 
complete equations of motion, the reduced non-linear equations of 
motion for the X-14 can be written as in Fig. A- 1-4. 

These equations can be put in state variable form by eliminating 
the derivative terms from the right hand side of equations (A- 1-10) and 
A-l-12). Rewrite the equations as 

I p = I r + L* 

X xz 

I r = I p + N* 

z xz 


(A- 1-1 3) 



where 


L* » (I - 1^) qr + I xz pq + T [-sin X ( z 2 cos o) 

+ Z 1 (v " V 1 V j ] + L P P + L 6 5 x + 

+ ^ V Sb (C t„ + % V 

a 

N* » (I x " I y ) pq ” I xz q r + 2 I e n £ q 

+ T [sin X (x 2 - sin o) - (v - ) /V ] + 

+ N r 6 + 1/2 p U 2 Sb (C + C 6 + C « ) 

6 z o n n 6 r n £ a 

z r a 

Upon solving Eqns. (A- 1-1 3) for p and r it is found that 


p = (I L* + I N*)/(I I 

r Z XZ X z 

r = (I L* + I N*)/(I I 

v XZ X X Z 


I 

XZ 


fc ) 



Figure A-l-4 
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.Summary of the Reduced Non-Linear Equations of Motion for the X-14 Aircraft 


u .« vr - wq - g sin 0 + - [ s ln a cos X - ~ (u - u ) ] 

j 


m 


•♦s'V 5 i< c i. 


p. off 


+ A C T ) sin ex - cos a] 


(A-l-7) 


v V P ■" ur y> cos 9 + — [sin X - ( v - v )] + pn 

m V v 2m U 0 Y 


(A-l- 8 ) f * 


w -« uq vp + g cos ^ cos 0 - ^ [cos a cos * + ~ ( w - w ) ] 

j 


w 


' is PU 0 2s [(C L 


(A- 1-9) 


(I - I ) 


p. off 
I 


+ ACL) cos a + C D sin a] 


P -= 


“T qr + ~I* ^ + Pq:) + T t" sin X (z 2 + £ - C0S c) + V ( v ~ V .J>] 


x 


L £ 

* i 2 p + r* «* + ii <= v Sb < c t + c , « > 

X X X *0 £ 0 


j 


w 


(A- 1-10) 


(I„ - I ) 

Z X 


, xz / 2 2. 1 t 

pr + y (r - ?)- - 21 n r + i- 

C e 1 


cos X 


sin o + cos a ) 4- ^ [x^ (w - w^) - (u — u^) ] 


M 6 


+ 1 *> + if 5 y + 2 T e u o Sc t ^ 


+ A C ) 

0 ft 

p. off 


(A-l-11) 


+ (c 


m 


+ * c m > a + fir* q + c_ 6j 


> off 


2U n 171 

0 q 


m 6 c 


(I - I ) 

* „x j. _jl* /• . 1 


f . yz ,, . l 

Pq + 7 -"* (p - qr) + ~2 I n q 
*_ lee 


N 


N 


+ j- [pinX (x - ; sin 0 ) - (v - v ) ] + q-^r + r~6 
z V j W l z l z 2 


] - 2 

+ - 57-0 U Sb (C + C 6 + C 6 ) 

0 n 0 n 4 r n 6 3 

r a 


(A-l-12) 


! 
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C. Numeri cal Evaluation of the X-14 Equations of Hot ion 

Numerical values for the constants in equations (A- 1-7) to (A-l-12) 

arc given in Table A- 1-2/ In addition, the following stability 

derivatives can be expressed as functional relationships [5]. 

2 

C. » 0.0030 cos 3/ deg 


-0.00014 cos Z 3/ deg 


C = 0.00118 cos 6/ deg 


n 


C ■ -0.0173/ deg 


HI, 


Furthermore, Figures (A- 1-5) through A-l-12) show graphs of C.^, 


C v , c » c m 
Y n~ m 


O' off’ Ver off’ \ 0Mcr o££ 


- V 1 C » 0 /T c’ 


AC /T and T vs. RPM. T is the non-dimensionalized thrust with T c 
L c 0 c 

T/ (1/2 U ‘ S). T is the thrust at the engine exhaust and T is the 
thrust at the end of the diverters. Hence, 

T “ T „ h 

where K a is the longitudinal diverter efficiency factor and ? x is the 
lateral side vane efficiency factor. Also, 

which demonstrate that the jet velocity is somewhat less at the exit 
of the diverters than at the exit of the engines. In order to reduce 
these graphs to a useable form, functional relationships must be ob- 
tained. It was attempted to fit all of the graphs with polynomial 
curves. This, however, was impractical for the more irregular func- 
tions and straight line approximations were used. The functions used 
are shown in Table A-l-3. 


J o 
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The reaction 

methods. The X-14 

«at Moffet Field. 

moment = >5. 

6 

y 


control nozzle effects wore evaluated 
's response to sinusoidal inputs was 
The second method is to note that 
<S = Force . Length 


by two independent 
analyzed by NASA 


or 



6 

y 

max 


where the length 


F 

max 
is the 


length 
dis tance 


to 


the discharge nozzle. 


F was 
max 


measured directly using transducers, and 5 was also measured di- 


y 

' max 


rectly. Numerical values associated with these calculations and a com- 
parison with the otner method are shown in Table A- 1-4. This agreement 

is quite good, especially when considering the errors and uncertainties 
involved in measuring quantities like the moments of inertia and the 
maximum force at the nozzle outlets. 


l 
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Tabic A- 1-2 

X-14 Constant Parameters 


P 

weight 

I 


2.378 x 10 _3 slugs/ft 3 


4182 lb. 


x 


2340 slug-ft 


= 3400 slug-ft' 


- 5400 slug-ft 
= 180 slug-ft 2 


mass center at station 99.84 in. and 0.58 in. right of center line 

x^ = 6 ft. 

x 2 = -0.0133 ft. 

z - 0.583 ft. 

z„ = 0.916 ft. 

2 

£. = 0. 167 ft. 


71 

e 

I 

e 

S 

c 

m 

b 


M /I 

q y 

l /i 

p x 

N /i 
r z 


m 


= 172S rad/sec at 100% RPM 
= 0.5 slug ft 2 
- 182.69 ft 2 
= 5.56 ft. 

«= 129.9 slugs 
= 33.83 ft. 


= -0.15 sec 
= -0.45 sec 
■* -0.20 sec 
- -11.4 


-1 


-1 


-1 


V 


-1 


» -0.0178 deg 
e 


. C = 0.576 rad 

"a 


M /I 
^ 7 

2 

= 0.0350 rad/scc /deg 

N /I 

Z 

2 

= 0.0166, rad/scc /deg 

L /X x 
^ X 

2 

= 0.0631 rad/sec /deg 

\ 

= 2000 ft/sec 


F = 0.91 

» 1.00 


Table A-l-3 


T 


T (Q) with H in percent of maximum RPM and T in pounds (for one 
°engine only) 


T = 44.15011 * ft - 1869.205 
o 

= 68.96552 * ft -A184.A83 
= 80 * ft - 5278. 

= 18.18182 * ft + 910. 


ft < 93.3 

i 

I , 

93.3 < ft < 99.1 
99.1 < ft < 100.1 
100.1 < ft < 101.2 


ft - n(T o ) 

ft = 0.02265 * T + A2.3375 
o 

= .0145 * T + 60.675 
o 

- .0125 * T + 65.975 
o 

= .055 * T - 50.05 
o 


T < 2250 
o 

2250 < T < 2650 
o 

2650 < T <2730 
o 

2730 < T < 2750 


k l /t c - f cy 


AC /T * -2.248398 x 10 ” 7 U 3 + 2.4S3009 x 10 -5 y - 1.332149 


L c 

x lo " 3 U 


0 


A V T C ' f ( V 


AC^/T = 4.900354 x 10 7 U Q 3 - 3.164620 x 10 J + 1.529505 


-5 „ 2 


x 10~ 3 U 


0 


S = W 


C D - .11 + l./(2. 152195 U Q + 12.5) 




0 


Si. 


<V 


""power off power off 

r *= -0.6 + 0.00237 U Q 

°power off = _ Q52 + 0>01A8 ( Uq _ 34 ) 


-0.27 + 0.001579 (U Q - 51) 


U 0 < 34 (ft/sec) 
34 < U Q < 51 
51 < U Q < 101 



M M 

“power off “power off 


<V 


M 


-4.1 + 0.04145 U 


“power off 


0 


- -2.7 + 0.1007 (U - 

*= - 1.0 


U Q < 34 (ft/sec) 


34) 34 < U Q < 51 

u 0 > 101 


c y « c y (8) 


C y = -1.8513 

= -0.35 + 0.393 (S 
= -0.26 - 0.362 (3 
- -0.55 


- 0.189) 

- 0.418) 


0 < g <0. 189 (rad) 
0.189 < 3 < 0,418 

0.418 < 3 < 1.22 

6 > 1.22 


C - C (3) 
n 0 n 0 

C = 0.113 8 


0 


0.022 

0.022 + 0.134S (g - 0.594) 
0.13 


0 < B < 0.1945 (rad) 
0.1945 < g < 0.594 
0.594 < g < 1.395 

8 > 1.395 


C S> = C £ (3) 

0 0 

C = -0.083 3 
*0 

=-0.014 + 0.0066 (S - 0.169) 
=-0.012 - 0.00865 (3 - 0.471) 
=- 0.02 

C L " C L 

power off power off (U^) 


L 


0.61 + 0.002763 U 


power off 


0 


0 < g <0. 169 (rad) 
0.169 < 3 < 0.471 
0.471 < B < 1. 395 
B > 1.395 


Uq < 51 (ft/sec) 


u o > 51 


0.75 
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D. Solution of the X-14 Equations of Motion for the Reference Control 
Variables 


The reference values for the elements of the control vector are 

desired for any given state vector x. Since a steady state reference 

condition is sought, set k *= (). Upon examination of the non-linear 
% 

equations of motion for the /X-14 (see Fig* A-l-4, page ), it is seen 
that u, v, and v are only functions of T, o and X . Therefore, an explicit 
solution for these variables should exist since there are three equations 
and three unknowns. To find this solution, first rewrite the equations 
as 


0 « A + T (sin o cos X - u ,/V, + f sin a) 

rel j 

0 = B + T (sin X ~ v ,/V.) 

rel J (A-l-14) 

0 «= C - T (cos a cos X + w . /V. + f cos a ) 

rel j 

where 

2 

A = in (vr - wq - g sin e) + 1/2 p U S (C sin a - C cos a ) 

U L/ r r JJ 

power off 

_ 2 

B = m (wp - ur + g sin $ cos 6) + 1/2 p S Cy 

2 

C = m (uq - vp + g cos $ cos 8) - 1/2 pU_ S (C 

1 power off 

cos a + Cp sin a) 


f " 4 C i/ T c 

After squaring, rearranging and adding equations, the result can be 


written as 

DT 2 + 2 ET - F = 0 
where 

D - 1 - 5 0 2/V J 2 - f2 + 2f < u rol 


- w , cos a) /V 
rel j 


sin a 
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E “ (A u , + B v . + C w ,)/V. - f (A sin a - C cos a) 
rel rel rcl j 

2 2 2 
F - A + B + C 

Hence, 

/ 2 

T = (-E + * E 3- DF )/D. 

Now since T is known, 

* E sin 1 (v . /V - B/T) 
rel j 

and now from Eqn. (A- 1-14) , 

cr c s i n ^ [ (u^^/V - f sin a - A/T)/cos ^ ] . 


(A- 1-15) 


(A-l-16) 


(A-l-17) 


(A- 1-18) 


Some additional information is needed to solve the moment equations 
since there are six unknowns and only three equations. Since both the 
three reaction controls and the three aerodynamic controls have approximately 
the same range of travel, let us assume for determining the control settings 
for reference flight conditions that the aerodynamic and the reaction con- 
trols are of the same magnitude, or that 

6 = 6 
x a 

6 =6 

y e 

6 = i 

2 r 

When the above numerically determined values of T, \ and a are 

substituted into the summary equations for the X-14, the equations can 

be solved for 6 , <5 ,6 , <$ ,6 , and 5 . The results are 
z e r x y z 

5 c 6 *= - { (I - I ) qr + I pq + T [-sin A (z. + l, cos o) 

xa y z xz nL 21 

+ z. v ,/Vl + L p + 1/2 p U n 2 Sb C. } /( L (A- 1-19) 

1 rel j p 0 o x 

+ 1/2 p U 2 Sb ) 

6 

a 
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6 - - { (I - I ) pr + I (r Z - p 2 ) - 2 I « r 

e z x r xz e e 


+ T [ cos X (z sin o+ x 2 cos o ) + (x j[ - Z;L u rcl >/V ] 


+ M q + 1/2 p U. Sc [ (C 

q 0 nip 


+ A C ) 


n v 


power off 


+ ,(C + A C ) a + - 

‘ m 2U m 

a ec a 0 q 

power off 


C q ] } / (M 


(A— ] 


+ 1/2 p IL Sc C ) 
0 m 6 


- - { (I - I ) pq - I qr + 2 I n q + N 
x y xz ^ e e r 


-- 2 


+ T [sin X (x 2 -t 2 sin o) - ^ v^/V ] + 1/2 0 U Q Sb (C fi 


+ C 6 ) } / (N . + 1/2 p U 2 Sb C ) 

n . a o U n x 

6 a 2 5 r 


(A-l 


-20) 


L— 2 1) 
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CHAPTER II 

LINEARIZATION OF NON-LINEAR EQUATIONS OF MOTION 
FOR A TIME INVARIANT DYNAMIC SYSTEM 


The non-linear equations of motion for a tine invariant dynamic 
system can be expressed as: 

X = _f ( x> u) (A-2-1) 

where x is the state vector, u_ is the control vector, and _f is a vector 
function of _x and u_ [1]. 

Suppose that for the system represented by Ean. (A-2-1) there is a 
given >c (t^ ) = end u_ (t) for t > t n , then the unique solution y (t, 


"O’ 


—0 


is determined and can be found. Now consider small perturbations c_x (t^) 
in and J"_u (t) in u (t). Hence, small perturbations in the solution 

$ 2Sq) can he expected and _x + c_x = _f (x + 6x f u. + 6 ij ) ; 6_x (t n ) - 6x 


0 — 0 

(A-2-2) 


Expanding the right-hand side, of Eqn. (A-2-2) in a Taylor series and 
retaining terms only of the first order, ar.d after substituting equation 
A-2-1 for x, ve have the variational equation 


<5x = A 6x + B cu; 6x (t n ) = ox. 

______ Q 


where 


A - 9f i 


b . . -!!i 

ij 3u 


(A-2-3) 

(A-2-4) 


with the partial derivatives being evaluated along the known solution, 
x (t) « ♦ (t, Xq) and control u (t). Hence, the notion of the dynamic 
system about a known path for a snail perturbation is seen to be governed 
by the linear ordinary differential equations (A-2-3) and (A-2-4), 
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The numerical evaluation of the partial derivatives in equation 
(A-2-4) can be accomplished by using the central difference approxi- 
mation. [2] A review of this method follows. The Taylor series for 

the function f. expanded about (x~ + Ax, u„) with 
x —0 — —0 


(Ax) T = (0 ... 0 


0 ... 0 ) 


is 


9f 


9f. . 

f i + A*. V = t ± (*o> i^) + ^ Ax. + 1/2 toj + ... (A- 2-5) 


a 2 f. 


with urg- ■ evaluated at (:<q, u_q) . Similarly, the function expanded about 
(Xq _ Ait, u Q ) is given by 


6f . 


3 2 £. 

1 


f i ( ^0 ' Z.' V = f i ( V V " x j + 1/2 ~2 2 ' X j ~ + ••• 

3 3x j 

3 

Since Ax is small, (Ax^.) and higher order terms are neglected. If 
equation (A-2-6) is subtracted from (A-2-5) , and the resulting equation 
solved for the partial derivative evaluated at the point (x^, u_ ) with 
Ax^ = 0 for k ^ j, the result is 


(A-2-6) 


9f . 

A. . = — ~ 


Vgo + AX. _V - h <Sp - **■ V 

2£x . 


(A-2-7) 


In a similar manner the elements of B are found to be 


B ij r 


8f i f i Ho + Au) - f. (x^, - Au) . 


9u . 


2(Au . ) 


(A-2-8) 


where (Au) = (0 0 0 ... Au. 0 . . . 0) . 

3 

A digital computer program that performs the above operations and 
its description can be found in Appendix A III. 
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CHAPTER III 

TRIM PARAMETERS AT3D STABILITY DERIVATIVES FOR THE X-14 

The numerical values for the trim parameters can be determined for 
any reference flight condition using equations (A-l-^4) through (A-l-21). 

I 

In addition, the stability derivatives can be calculated using the lin- 
earization method of Chapter II and the equations of motion for the 
X-14 of Chapter I. 

Trim parameters and stability derivatives for thirty flight conditions 
are shown in Tables A- 3-1 to A- 3-30. These were calculated on a digital 
computer using the FORTRAN IV programs of Appendix V. The reference 
conditions are for wings level and forward velocities ranging from -10 
feet per second to 4-169 feet per second. The verticle velocity ranges 
from -10 to 4-10 feet per second and the horizontal lateral velocity goes 
from zero to 415 feet per second. 

The zero forward velocity condition is represented by a numerical 
value of 0.01 feet per second which was necessitated because of the non- 
dimensionalization method used. If u = 0, then certain segments of the 
equations of motion could require division by zero. Any error that might 
be introduced by this approximation is trivial for the conditions being 
investigated in this paper. 

The "A Matrix” represents the relationship between the tine derivatives 
of the state variables and the state variables as listed in the reference 
condition column. The "B Matrix” is the relationship between the state 
variables and the control variables as listed in the trim parameter column. 
The "Wind Disturbance Matrix" represents the effect of wind disturbances on 
the state variables. The columns of the wind matrix are for the X, Y and 


Z directions, respectively. 
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